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Abstract 

<N ■ 

^ ' We analyze a thought neutrino oscillation experiment in which a beam of neutri- 

nos is produced by electrons colliding with atomic nuclei of a target. The neutrinos 
I . are detected by observing charged leptons, which are produced by neutrinos colliding 

I with nuclei of the detector. We consider the case when both the target and detector 

■ nuclei have finite masses. (The case of infinitely heavy nuclei was considered in the 

. literature earlier.) 

p i' 

' : 1 Introduction 

(D . 

1^ • Despite an impressive number of theoretical papers published during the last 40-50 years, 

>• ■ the phenomenology and description of neutrino oscillations is still a subject of heated 

^ . debates. In particular, there is no consensus on the assumptions of equal energies or 

equal momenta of the three neutrino mass eigenstates uj, j = 1, 2, 3. 

The equal momenta scenario was introduced in a pioneering paper on neutrino oscilla- 
tions by Gribov and Pontecorvo PP, used by Fritzsch and Minkowski [2] and then by many 
other authors. The equal energies scenario was presented by Kobzarev et al. jH], who 
considered all three virtual neutrinos produced by a monochromatic beam of electrons 
colliding with infinitely heavy nuclei (M — > 00). Since the recoil energy of such nuclei 
is zero, all three neutrinos have equal energies, the same as the energy of the electron. 
Stodolsky Lipkin [3] and Vysotsky presented general arguments in favor of equal 

*e-niail: dolgov@itep.ru 
''^e-mail: okun@itep.ru 
•l-e-niail: nirotaev@niail.ru 
^e-mail: schepkin@itep.ru 



1 



energy scenario for realistic thought experiments. Still, in the most recent and authorita- 
tive review of particle physics in the contribution by Kaiser [Tj neutrinos oscillations are 
discussed on the basis of equal momentum scenario. Note that the same attitude one can 
find in his previous review |HI, while in 2000 |^ both equal energy and equal momentum 
scenarios were considered on the same footing (all this - in the oversimplified "neutrino 
plane wave approximation"). 

In refs. fT^ , [2^ , [7J , JHj , [D] plane wave free neutrinos traveling from the production point 
A to the detection point B were considered without discussing their progenitors. We 
will refer to such descriptions as reduced ones. Following the argument of ref. jH] it is 
evident that in the "reduced approach" only the equal energy scenario is self-consistent. 
Otherwise the neutrinos are produced at the point A not in a given flavor state, but in a 
state whose flavor oscillates with time. 

In ref. jS] the progenitor (electron) was described by a plane wave. There exists a vast 
literature in which both the progenitor and offspring particles are described not by plane 
waves but by wave packets (see review by Beuthe fD])- In the present paper our attention 
is concentrated on the incoming and outgoing particles. We will find that depending on 
the properties of these external particles both differences of energies and momenta of 
different neutrino mass eigenstates are non-vanishing, but the energy difference is much 
smaller than the momentum difference when the energy transfer to the target nucleus is 
small. 

In this note we are going to consider a more realistic situation than in ref. |3] , namely, 
when the beam of electrons is not monochromatic (it is described by a finite-size wave 
packet), and the mass of the target nucleus is finite. Now the recoil energy of the nucleus 
cannot be neglected. The detection of neutrino occurs when it interacts with another 
nucleus (in detector). When a nucleus is in a crystal it is described by a wave function 
with a characteristic momentum spread about 1 keV and vanishing mean momentum. 
When the nucleus is in gas, its momentum is not vanishing, while momentum spread 
is smaller. (Note that even for an infinitely heavy nucleus in crystal the spread of the 
momentum is non-zero). 

We will prove that in the case of finite masses of nuclei A or B neutrino oscillations 
disappear in the limit of the vanishing momentum spread of the electron wave packet 
(plane wave limit). 

The structure of this paper is the following. We introduce "little donkey" diagram 
with a virtual neutrino propagating between production and detection points in section 2. 
The amplitude is derived and analyzed in section 3. The expression of the phase difference 
responsible for oscillations is discussed in section 4. The probability of neutrino oscillations 
and their suppression is discussed in section 5. Section 6 is devoted to concluding remarks. 

2 Probability and "little donkey" diagram 

Let us consider interaction of an electron e with a nucleus A of mass Ma in a target. 
Neutrino produced in this interaction collides later on in a detector with a nucleus B 
of mass Mb and produces a charged lepton /. As a result the whole process looks like 
e + A + B ^l + C + D, (see Fig. 1). 
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Figure 1: Little donkey diagram 



The electronic neutrino v^. produced on nucleus A is a superposition of three neutrino 
mass eigenstates: Vf, = Ueit^i where Ui is the state with mass mj. Each mass eigenstate 
propagates independently between nuclei A and B. Interaction with B results in projec- 
tion of the three neutrino propagating states on the state z/; = Ylii^n^i- Here U is the 
unitary mixing matrix, the first and second indices of which denote respectively flavor 
and mass eigenstates. 

In oscillation experiments the nuclei C and D are not registered, while the energy and 
momentum of the lepton / are measured with low precision. Thus the probability of the 
whole process is obtained by integration of the amplitude squared, possibly weighted with 
the detector resolution function, over dTpidTpcdpo 



P{Pe,PA,B, Ma,B,^A,B, 0-A,b) = 



X Pij{Pe,l, Pa,b, Ma,b, ^a,b, (Ta,b, Pc,d, ^^'^)^^^^^- 

Here Pij = A*Aj, where Ai is the amplitude for a given neutrino state with mass rrii] 
symbols like Pa,b denote p^i, p^. Let us express Aj in terms of the wave functions of all 
interacting particles: e. A, B, C, D, I and neutrino propagator G{xi,X2) : 

Aj(pe,i,PA,B, Ma,b,^a,b,(^a,b,Pc,d, Ec,d) = J rfxidxa^tidta X 

X ?/'e(xi, ti)^A(xi, tl)V^B(x2, (xi , t, X2, t2)V^r(x2, t2)^C (^1 ' ^1 )^Z) (^2 , ^2)- (2) 

In ref. ^ electron e and lepton / were described by plane waves: 

while nuclei A and B were described by (5-functions in configuration space ^. 

"'^The equality signs in equations throughout the paper should be taken "with a grain of salt" because 
we omit some obvious normalization factors. This makes the formulas easier to read without influencing 
the physical results, e.g. the ratio of oscillating terms to the non-oscillating ones. 
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Now we assume that nuclei A and B are described by the finite-size wave packets: 



where and are the central momenta of the wave packets of nuclei A and B respec- 
tively, while qyi and qs are the corresponding running momenta. In eqs. (jH) and below 
we use Gaussian wave packets. The main features of our results do not depend upon the 
specific form of the packets. In our subsequent publication we are going to take for them 
a general form. 

The wave functions of nuclei satisfy Klein-Gordon equation as we consistently neglect 
the spins of all particles. All external particles are assumed to be free and hence their 
momenta satisfy the on-mass-shell condition: 



EM = V^l + Ml , (5) 

where the index n denotes A, B and e. The same is true also for outgoing particles. 

Since in neutrino oscillation experiments the final nuclei C and D are not registered, 
one can choose for their wave functions any basis in the Hilbert space, and then integrate 
the probability over all Hilbert space. In what follows we take plane waves as the complete 
and orthogonal set of wave functions: 

^c(xi,ti) = e^Pcxi-*i?c(ti-tA) ^ 

where E"^ = p'^ + M^ and Ejy = p|, + M|,. We think that description of outgoing particles 
by wave packets in the amplitude is not a consistent procedure, because generally wave 
packets do not have orthogonality property, neither they form a complete set of functions. 

Here we would like to touch upon a subtle point. In a more or less realistic thought 
experiment the target and the detector are solids, therefore the on-mass-shell condition for 
nuclei is only an approximation. This approximation seems to be reasonable for ordinary 
matter, where nuclei are weakly bound. 

Instead of the plane wave of ref. the wave function of the electron is described 
now by the one dimensional wave packet with definite direction e = Pe/|pe| of the beam: 

M^uti) = j e"^'''^e^^=-(^^-^^)-*^=('i^)(*i-*=)ciqe. (7) 
Here and in the following: 

d(\e = T-^WedqedVLe, (8) 

iQe I 

where Vie is the corresponding solid angle, and Qe = |qe|- We will show below that the 
oscillation terms vanish when ae tends to zero. We choose the one-dimensional packet 
only because of technical simplicity. The result can be obtained in a more general case. 
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3 Neutrino Green function and the amplitude 



Following ref. [H] we replace the neutrino Green function with the propagator of a scalar 
particle of mass rrij, where j numerates neutrino mass eigenstates, j = 1,2,3; it is clear 
that fermionic nature of the neutrino (as well as of e and /) is not essential in the problem. 
Thus 



G',(x,t) = - — - / — — — . dkduj . (9) 



For each za- the amplitude of the process is written as 



/ (pe-qe)^ _(qA-PA)^ 
g 2^1 e^^^^^'^~^^'^~^^^^^'^~*-<='^dc[ee ^ g«qA{xi-XA)-i-BA(qA){ti-tA)^q^ 



X 



X 



^ g-ipcxi+*Bc(il-tA)g-*PDX2+iBB(t2-tfl)^XirfX2rftirft2 = 

(Pe-qe)^ _(qA-PA)^ _(qfl-PB)^ 

2^ giqe(xi-Xe)+j-E;e(te-iA)g 2^7^ gjqA(xi-XA)g 2^ ^ 



X g«qs(x2-Xs) ]_ g-«fcjlxi-X2l+iaJ-(tA-ts)g-jPiX2g-ipcXl-jpDX2^-j^^^j^^^ 

47r|xi - X2I 

X S{E, + E^iqA) + EsiqB) - Ec - Ei - ED)dqJqAdciB , (10) 
where we use dtidt2 = ^d(ti + t2)d{ti — t2), and 

! e-^-(*^-*^)G,-(Xi - X2,ti - t2)d{ti - t2) = ——^ ^g.v/^^|x,-X2| 

J 47r|xi-X2| 
The parameter uj is defined by 

u = uj{qA,qB) = E, + EA{qA)-Ec = Ei + Ed - EeiqE), (12) 

and 



k^^^u^-m]. (13) 

Though kj looks like a three-momentum, in fact, it is a short-hand notation, usually 
arising in description of propagation of spherical waves with definite energy. 

The integration over d{ti + 12) in eq. ()10|) gives 5-function leading to energy conserva- 
tion. 

The further analysis of the problem is greatly simplified if the distance |x^ — x^l is 
much larger than the sizes of wave packets of nuclei A and B. To take this into account, 
let us shift the variables of integration: 



Xi = Xa + X' 



X2 = Xb + X2 . 



(14) 



5 



The wave packets of the nuclei A and B are essentially different from zero if |x'^|, 
Ixgl <^ |xa — x^l, hence 

|xi - X2I ~ |xa - xsl + (x'l - X2)n, (15) 

where n is the unit vector in the direction Xyi — x^. 

By substituting eqs. (fT^ and (|THjl into eq. (fTUjl we obtain after integration over dx'^ 
and dx.2: 



\n=e,A,B J 



~Pc)xA-«(Pi+Po)xs-«qeXe+jBe(qe)(te-tA) y- 



(Pn-qn) 

X e ^"--•^•^ ^^-A 6{qe + qA- kjU - pc)S{qB + k^n - pi - p^) x 
^ exp {zfc,|xA - x^l + .^,(tA - t..)}^q^^q^^q^ ^ 



E ikj |x^ -xs |+ja;(qAj ,qBj )-{tA-tB) 
En\qn',) + ^e(ge) - " " ^d) ^ ^ X 



_ (P„/-q„/j)^ (Pe-qe)" 

X e*(Ps"Pc)xA-i(P!+Pi3)xs-ipeXe + i_Be(ie-iA)g ^n' = A,S ^^^^ 



The j-dependent momenta q^j and q^j are defined as: 

qAj = kjxv + PC - gej-e , 
qi?i = Pi + Vd- kjU . 

We integrate over dqe using (jSJ. The result is: 



(17) 



dqe J 47r|xA - x^ 

X 6*^^"=^ 



(pn-q„,) 



i(lJeje-pc)xA-i(p!+PD)xs-iqeXe+i-Ee(qe)(ie-tA)p En = e,A,S 
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where Uj = uj{qAj,qBj) (see eq. ^s(ge) = Y.n'=A,B ^n'{qn'j) + Eeiqe)-Ec-Ei~ED, 

and {dEj:,{qe) / dqe)~^ is the Jacobian, left after integration of the energy 5-function in (fT^ . 
(In eq. p6|) there are three (5-functions, one of them expressing the energy conservation, 
while the other two refer to momentum conservation in A and B vertices.) 

We have already stressed that kj is not a momentum, but a parameter characterizing 
spherical neutrino wave. Now we see that in the case of very large distance |xa — xb| 
the parameter kjU does play the role of the neutrino momentum. We are faced with the 
situation when neutrino being virtual particle at short distance from the source becomes 
effectively real at large distance, near detector. 
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4 Phases of the amphtudes 



We are interested first of all in the phases of amplitudes of the process considered. From 
eq. (fTR|) one can see that the phase of Aj equals to 



= kj\yiA - xb| + ujj{tA - ts) + {qeje - Pc)xa- 

- (P/ + PZ))xb - gej-eXe + Ee{Pej){te " ^a), (19) 

and dependence of kj, Qej and ujj on rrij is given by the system of equations (fT^ . (fT!?|l . 
()17|). and by the on- mass-shell conditions for nuclei and electron. 
From eq. (fT^ it follows that 

4>ij = 4>i- 4>j = |xb - ^A\{ki - kj) + {Ui - UJj)(tA - tB) + 

+ {qei - qej){^A - Xe)e + {Ee{qei) - Ee{qej)){te - tyl)- (20) 

In equation (pUj) the difference of the electron energies could be expressed through the 
difference of the corresponding momenta: 

dE (q ) 

[Eeiqei) - Eeiqej)) = -^^{qei - qej). (21) 

dqe 

Since neutrino masses are much smaller than energies and momenta of external parti- 
cles^ we may write: 



dk ■ 

5kii = kj — /co ~ imj — m^^ •' 



duj 



2uo ■ (1 - vb ■ n 



2 2 

^ (22) 



mj=0 



5uJij = uJi- ujj ~ {m,- - m ■) ■ 



3/ B ^23) 



2 2\ ""^ej 



gej ^ (m^ -m .) 



(i(m' 



„^|=o 2a;o 1 - VB ■ n 

(mf - m|) (vb - VA)n 



(24) 



where v„ = q„o/-E'n(qno)) and the subscript denotes e, A, 5. 

The quantities q^o, -E'n(qno) and cuq are defined by the external parameters from eqs. 
(fT^ . (jlSI), and (fTTjl at = 0. In particular: 

(pD-pO^ + M|-(Eg + EO^ ^ (PC - Pe? + Ml - (Ec - i^e(Pe))^ .25) 

""^ 2[(pz) + pz)n - Ed - ^i] 2[(pe-pc)n + Ec-Ee(Pe)] ■ ^ ^ 



^Let us point out that in the hmit Ma,b ^ co the neutrino energy tends to: uji) + Ea — 

Ec = El — Eb + Ed = E^ + 0{1/Ma,b), the value defined by the energy conservation in the process 
e + A + B^l + C + D, and since (qn)/_E(q) 0, the phase difference approaches its standard value 

? ? 9 9 

— rri: — m, 

-|xi3 - xaI = 2E'^^^ ~ 
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In eq. I^E^ the sign "=" means exact but somewhat useless equahty, because p^^ and 
Ed are not measured, while p/ is measured with low accuracy. As for the sign it will 
be used in what follows because Pe and Eg are known and essentially define the value of 
ujo: 

uJo = E, + 0{^). (26) 
From eqs. pi^ . it follows that 

2 2 

— ml. 

4>ii = ; — rX 

2cuo(l - VB ■ n) 

-|xb - x^l + {tB - tA)^B ■ n + [(xe - ^a)g - \w,\(t, - Ia)] ' ^ • (27) 

It is convenient to choose the parameters tg, t^, Xg, and x^ in such a way that Xg = x^ 
when tg = ^A- This convention corresponds to the classical picture of eA-coUision and 
allows to simplify eq. ()27|): 

2 2 

m| — m- 

^ij = n — 71 — V i-\^B - xaI + {ts - tA)^B " n) . (28) 

ZLLJo[i — vb ■ nj 

From eqs.(j2Hl), (EH) and ^ it follows: 

4>ij = 5/Cij|xB - Xa| - 6uJij{tB -tA) = (29) 



2 2 

mf — m 



-|xb - xa| + 6uJij [|xb - x^I - (ts - Ia)] . 



2ujQ 

The first term in the phase is the standard phase of oscillation theory, while the second 
one is an additional term which depends upon the size of the electron wave packet. 



5 Neutrino oscillations and their suppression 

By using eq. (fTUI) we obtain the following expression for describing neutrino oscillations 
in the right-hand side of eq. ((H): 







'A 


2 




2 






- x^l 


2 



dqe ) 

xexpf- ^ ^5^^^- E ^^^^] . (30) 

\ n=e,A,B " n=e,A,B " 



where (pij is given by eq. 

This formula allows to compare the oscillating terms {i 7^ j) with non-oscillating {i = 
j) ones, and thus to analyze the strength of oscillations as a function of the momentum 
spread of the electron wave packet. 
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For easier comparison with ref. [3] we assume in what follows that is much smaller 
than a A and as- Let us define 



h - ^^P(- 2^^' ^ ^^^^ 

and assume that |pe — qei| < |Pe — qe2|, |Pe — Qesl, then /i ^ /2, /s in the limit of vanishing 
CTg. The leading diagonal term 

~/i'>P22,P33. (32) 

Comparing Pn with the non-diagonal terms we conclude: 

Pii~/i'>Pi2~/i-/2, and Pn>Pi3,P23. (33) 
Considering the ratio 



\2 



exp I ■ (qe, - qe.) 1 (34) 



fj I Vf e 2 

= exp 



(35) 



and using eq. one finds the crucial parameter of suppression to be 

Qej - Qei _ (vj? - va) • n 

where 

is ij oscillation length. 

The Gaussian factor in eq. (jnn|) makes it obvious that for |q„i— qrij| ^ cr„, the oscillating 
terms become exponentially suppressed in comparison to the non-oscillating ones. 

In conclusion of this section let us make the following remark. Though the suppression 
for vanishing ae is obvious, it is clear, that in "a realistic thought experiment" -C Lose, 
and hence the suppression is very weak. 



6 Concluding remarks 

1) We see that the alternative "equal energies versus equal momenta" is naturally resolved 
if one consistently uses the standard rules of quantum mechanics and in particular quan- 
tum field theory. In the example, which we consider here, using the propagator of virtual 
neutrinos and mixed description of initial (wave packets) and final (plane waves) particles 
all kinematical variables are uniquely defined. In particular, when we go beyond plane 
wave approximation for initial particles there is no equal momenta nor equal energies. 
However, still \ujj — uJi\ -C \kj — ki\ at least for non-relativistic nuclei. Similar conclusions 
were obtained in refs. ^21, ^B] for oscillating neutrinos produced in pion decay (see also 

2) For the plane wave of the initial electron and finite mass nuclei, neutrino oscillations 
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disappear unlike the case of infinitely heavy nuclei. For a finite but small momentum 
spread of the electron wave packet, the neutrino oscillations are suppressed. 

3) For realistic parameters of the electron wave packet the above suppression is small and 
therefore can be disregarded. 

4) The Green function used to describe neutrino leads us to the situation when the neu- 
trino being a virtual particle at short distances from the source, becomes effectively a 
real particle at large distances, near detector. This is the standard case in the scattering 
theory. 

5) With localized "meeting points" eA and IB the time dependence of the oscillation 
probability is not essential. (The time moments and ts enter the expression for (pij 
with small coefficients proportional to velocities of nonrelativistic nuclei.) 
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